Searches on star graphs and equivalent oracle problems 
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We examine a search on a graph among a number of different kinds of objects (vertices) one of 
which we want to find, fn a standard graph search, all of the vertices are the same, except for 
one, the marked vertex, and that is the one we wish to find. We examine the case in which the 
unmarked vertices can be of different types, so that the background against which the search is 
done is not uniform. We find that the search can still be successful, but the probability of success 
is lower than in the uniform background case, and that probability decreases with the number of 
types of unmarked vertices. We also show how the graph searches can be rephrased as equivalent 
oracle problems. 

PACS numbers: 



I. INTRODUCTION 

One of the main reasons quantum walks were originally 
introduced into quantum information was the hope that 
they might provide a useful conceptual framework for the 
development of new quantum algorithms [IJ . This hope 
has been fulfilled [2[-[5|- In addition, they have recently 
attracted considerable experimental interest. They have 
been realized in optical lattices @, 0] and with photons in 
waveguide lattices 0. Analogs of quantum walks, which 
are based on the wave nature of classical light, make use 
of the fact that the most important element in the behav- 
ior of a quantum walk is interference. These were 
posed in [| E3 

and experimentally realized in Ref. 
The use of interferometers to realize quantum walks was 
proposed and recently experimentally implemented 
[12( | . Quantum walks with trapped ions have been theo- 
retically described [HI, [l4| and experimentally realized as 
weii \m [l6j . Finally, experiments in [l7| and [l8j provide 
a framework in which walks with more than one walker 
can be studied. 

Here we will look at quantum walks on star graphs 
and equivalent oracle problems. A star graph has a sin- 
gle central vertex with N edges emanating from it. Each 
of these edges is connected to a single external vertex. 
Thus the graph resembles spokes extending from a hub. 
We shall denote the central vertex by and the exter- 
nal vertices by 1 through N. Star graphs form a natural 
arena in which to study quantum searches. Some of the 
external vertices may have different properties from the 
others, and our task is to find these vertices. In a stan- 
dard quantum search one object has different properties 
than all of the others, all of which are the same, and 
which constitue a kind of background. Here we want 
to examine searches on more complicated backgrounds. 
Suppose that we want to find an object with a particu- 
lar property, but the other objects that we must search 
among are not all the same. In the case of a star graph, 
this could be realized by having the external vertices re- 
flect the walking particle with one of three phases. One 
vertex (the marked vertex) reflects with the phase we 
want to find, and each of the other vertices reflects with 



one of the other two phases. As we shall see, it is possi- 
ble to find the marked vertex using a quantum walks, but 
there is a limit to the probability of successfully doing so. 

We shall employ the scattering model of a quantum 
walk in which the particle resides on the edges of the 
graph and at each time step scatters at the vertices [lj| 
l2fj |. Each edge has two states, corresponding to the two 
directions a particle on it can be going. In the case of a 
star graph, the edge between vertex and j is associated 
with the state |0,j), which corresponds to the particle 
being on that edge and going from to j, and \j, 0), 
which corresponds to the particle being on that edge but 
going from j to 0. The states corresponding to the edges 
form an orthonormal basis for the Hilbert space of the 
particle on the graph. 

In addition to a Hilbert space we need a unitary oper- 
ator that advances the walk one time step. We shall give 
a detailed description of this operator shortly, however 
when a particle scatters at an external vertex the effect 
of the operator is to reflect the particle with a phase, 
and the phases can be different at different vertices. It 
is these phases that determine the behavior of the walk. 
One case is where all of the phases are except for one 
which is 7T. This simply reproduces the standard Grover 
search, and the particle becomes localized on the edge 
with the 7T phase shift after 0{\/N) steps. As we men- 
tioned, we can do a search on a more varied background 
by allowing the phases to have more values. If roughly 
half have phase shifts of 2tt/3, roughly another half have 
phase shifts of — 2ir/3, and a few have phase shifts of zero, 
the particle will become localized on the edges with phase 
shift after 0(^/N) steps, but the maximum probability 
of finding it there is not close to one, as in the case of the 
Grover search, but 3/4. Adding further variety to the 
phases causes this probability to decrease. 

We will begin by examining the eigenvalues of the uni- 
tary operator that advances the walk one step. One result 
is that if most of the phase shifts are the same and only 
a few are different, the eigenstates group themselves into 
two subspaces, a big, relatively boring, subspace in which 
not much happens, and a small space in which interest- 
ing things can happen. This allows one to focus on what 
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is happening in the small space and to ignore the large 
one, which has the effect of simplifying the analysis of 
the walk. 

Next, in order to connect the quantum walk descrip- 
tion of a search to the usual oracle-based one, we analyze 
what happens when the function the oracle evaluates is 
not a binary- valued one. We examine two different uses 
of the oracle. The first way of using the oracle reduces the 
problem to a standard Grover search, and the probability 
of finding the marked vertex is close to one. This use of 
the oracle, however, is not equivalent to what happens in 
our quantum walk. We, therefore, explore a second way 
of using the oracle, which is equivalent to the quantum 
walk on the star graph. Using the oracle in this way, for a 
function that can take more than two values, corresponds 
to allowing more than two phases in the quantum walk. 
For this use of the oracle, we further study how increasing 
the number of possible function values affects the num- 
ber of steps in the search and the maximum probability of 
finding the marked element. These results carry over im- 
mediately to the quantum walk search on the star graph 
with multiple phases and give us the number of steps re- 
quired to localize the particle and the extent to which it 
can become localized. 



II. EIGENSTATES AND EIGENVALUES OF 
STAR GRAPHS 

Let us now specify the details of the walk on a star 
graph. As noted in the introduction, if the star graph 
has N edges, we denote the central vertex by and the 
other vertices by 1 through N. When the particle making 
the walk is reflected from the vertex at j, it picks up a 
phase factor e lc ^ j . The action of the unitary operator that 
advances the walk one step is given by 



U\0,j}=e^\j,0), 



N 



U\j,0) = -r\0,j)+t l°' fc > 

N 

= -io,j> +*!>,*>, (i) 



fc=i 



where r = (N - 2)/N and t = 2/N, Let 



N 



3=1 



(2) 



If \ip) is an eigenvector of U with eigenvalue e z9 , i.e. 
U\ip) = e l6 \il>}, then we have that 



JV 



fc=i 



(3) 



Solving the second equation for j3j and inserting it into 
the first gives 



tS 



(4) 



where S = J2f=i a j- There are now two cases. First, 
suppose S 7^ 0. Then solving the above equation for ay 
and summing it over j gives 



JV 



1 



3=1 v ' 

so that the eigenvalues are determined by the equation 

1 



(5) 



1 = t 



N 

E 

3 



1^ ( e i(2fl-^) + ]_) ■ 



If S = we have that 



which implies, if otj ^ 0, that 



(6) 



(7) 



(8) 



The only way this can be satisfied is if two or more of 
the d>i's have the same value. 



Let's look at an example. Suppose 



7r and 







for j > 2. This is, of course, the standard Grover case. 
Setting z — e 2%9 Eq. (J6]) becomes 



2rz + 1 = 0, 



(9) 



giving z — r±i(l — r 2 ) 1 / 2 . Setting e 2l6a = r + i(l — r 2 ) 1 / 2 , 
we find that this gives a value for 9q that is of order 
1/VN. We thus have four eigenvalues of this type, ±e l6, ° 
and ±e~ l6 °. Eq. © yields the eigenvalues ±i. Each of 
these eigenvalues is (N ~ 2)-fold degenerate. It is also 
quite straightforward to find the eigenvectors. For e l0 ° 
we find 



\ti)=ri 



1 



1 -e 



-(|l,0)-e' e °|0,l» 



•ET-^(i^ ) + eieo i ^-)) 

3=2 + 



(10) 



and for — e i6 ° 
16) = V 



1 



I _ e 2i6 
N 



(|l,0)+e^|0,l)) 



3=2 



(11) 



The expression for the eigenvector (£3), which corre- 
sponds to eigenvalue e~ l0 °, is the same as that for 
except with 9 Q replaced by — # , and the expression for 
I £4), which corresponds to eigenvalue — e~ l6 °, is the same 
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as that for |^) except with 9q replaced by — 6$. The 
normalization constant 77 is given by 



7] = 



sin(20 o 



1 



[N- (N -2)cos(28 )} 1 / 2 y]v' 



(12) 



by \£,i) an d I £3)- After m steps the state of the particle 
is 



\4, m ) = _( e ^0| 6)+e -^ | 6))+O(7V -l/2 ); (15) 



The eigenvectors corresponding to the eigenvalues ±i are and when m8o = ir/2, we have 



just determined by the condition ~^2f =2 a j — 0, so there 
is a great deal of freedom in choosing them. One possible 
choice for the eigenvectors corresponding to +i is 



IC 



W) = V2 ^_ 1} E ^-^/("-D^ 0) + i| 0) ,"», 

for m — 1,2, ■ ■ ■ N — 2. The eigenvectors corresponding 
to —i can be taken to be of similar form. If we start the 
walk in the state 



IV>o) 



1 



>2N 



N 

E 



(L?,0) + |0,j», 



(14) 



we find that |V> ) = + |6))/V2, assuming iV > 1, 
so that the whole walk takes place in the space spanned 



l^)--^=(|l,0)-|0,l)) + O(^V- 1/2 ). (16) 

Therefore, to good approximation, the particle is located 
on the edge connected to the vertex with the tt phase 
shift, just as in the Grover algorithm. 

Now let us move on to a more complicated situation. 
Suppose that each of the tj>j 's takes on one of three values 
4>\, <f>2 or 03, and, in particular, n\ are <f>±, n% are 4>2, and 
nz are 03, where m + n 2 + n 3 = N. Then Eq. ^ gives 
us 



N 



2 ^ ze-^i + 1 ' 
3=1 



E 



(17) 



which becomes 



E; 



z 2 + 2 



(18) 



where xj — rij/N. 

One case in which the equation can be solved is if <pi — 
2ir/3, 4>2 = —2ir/3 and </>3 = 0. In addition, we shall 
assume that N is odd, and n\ = n 2 = (N — l)/2, and 
ri3 = 1. In that case, the above equation becomes 

The solutions to this equation are z — 1 and 

z = - 1 + 4 ±i (|-i^) ,/2 - < 20 » 

Note that the last pair of solutions for z yields eigenvalues 
of the form ±[i ± 0(N^ 1 ^ 2 )]. Comparing to the previous 
case we can surmise that this situation will again lead to 
a search in which we can find the vertex corresponding 
to the phase of zero in 0(y/W) steps. This is, in fact 
the case, but rather than use the exact eigenvalues to 
show this, we shall make use of an approximation method 
which allows us to find a solution in a more general case. 

We want to look at the situation in which one of the 
Xj's is small, typically of order 1/N, and N 3> 1. In that 
case the above equation can be solved perturbatively. If 



the full third degree polynomial above is f(z), then let 
fo(z) be the same polynomial but with the terms of order 
1/N set equal to zero, that is, fo(z) is the N — > 00 limit 
of f(z). We then get a zeroth order solution to f(z) = 0, 
which we shall denote by zq, by solving fo(z) = 0. 

We can then get corrections to the zeroth order solu- 
tion, by setting z = zq + 5z substituting back into the 
cubic equation, keeping the lowest order surviving terms, 
and then solving for Sz. This will give us 

= f(z Q + Sz) 
-f(z Q ) + [f[ ) (z ) + O(l/N)]5z 

+ [f^z Q ) + 0(l/N)](Szf + --- . (21) 

Note that we have substituted fo + 0(1/N) for / in the 
derivative terms, e.g. f' {z Q ) + O{l/N) for f'{z Q ), because 
the fo terms are typically of order one. Because f(zo) is 
usually of order 1/N and the (Sz) 2 term can be neglected 
in comparison to the Sz term, we will usually have Sz ~ 
1/N. This situation is not interesting, and let us see why. 
If z = e i(>0 , and z + Sz = e ^ e ° +se \ then Sz - 1/N will 
give us 80 ~ 1/N. After m steps we will have a factor 
of e lm ( e o+<50) multiplying the eigenvector. The part of 
this phase factor that contains the information about the 
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unusual vertices (the ones with the rare phase shift) is and 
e imS9 anc j ^ n j s phagg f ac tor will differ substantially from 
one when m is of order TV. That means that in this case, 
the unusual vertices would have a significant effect on 
the dynamics after N steps. However, we could find the 
unusual vertices classically just by checking each vertex 
to see if its phase shift was one of the rare ones, and this 
procedure would take iV steps too, so we have gained 
nothing. 

What we really want, to get some kind of a quadratic 
speedup, is for SO ~ 1/y/N. Examining Eq. (|2ip . we 
see that this can happen if the order one part of the 
coefficient of Sz, which is /q(z ), vanishes. We then get 
a quadratic equation for Sz 



= f(z Q ) + [0(1/N)]5z + f^'(z )(5z) 2 



(22) 



which typically results in Sz being of order N~ x / 2 . The 
condition /q(zq) — implies that zq is a double root of 
the equation fo(zo) — 0. Summarizing, in order to find 
a situation in which we can have a quadratic speed up, 
we want the solution of the zeroth order problem to be a 
double root. 

Let's give an example where that happens. We are 
going to assume that almost half of the phases are <fii = 
27r/3, almost half are <j>2 = — 27r/3 and the remainder are 
<f>3 = 0. We will take most of the phases from j = 1 to 
j = N/2 to be 27r/3, and most of the phases between 
j = {N/2) + I to j = N to be -2tt/3, and we will split 
the phases that are zero putting some between j = I to 
j = N/2 and some between j = (N/2) + 1 to j = N. 
In particular, let us assume that for j — 1, • • • ,71.31 and 
for j = {N/2) + 1, • ■ • , (N/2) + n 32 the phases are 0, for 
j = "-31 + 1, • • • , (N/2) the phases are 27r/3 and for j — 
(N/2)+n 32 + !,-■■ , N the phases arc -2tt/3. Therefore, 
the total number of edges with phase is n 3 — ti 31 + n 32 , 
the total number with phase 27r/3 is n\ = (N/2)— 7131 and 
the total number with phase — 2-7r/3 is n 2 = (N/2) — 7132. 
We will consider the case where ^3 <C N. 

Let us first get the zeroth order solution to Eq. (TT51) . 
To do so, we neglect all quantities of order n%/N , which 
means we set x\ = x 2 = 1/2 and x 3 — 0. Doing so, we 
get the equation 



Xie i(4>2+4>3) + X2e i (<t>i+4>: i + ^e^i + fe) 
-2tti/3 



= ~2 + *3i(l - + 132(1 - e 2 " /d ). (25) 

Upon making our substitution and keeping only quanti- 
ties of up to second order, we first find that, as expected, 
the term proportional to Sz vanishes, and we are left with 

(Sz) 2 - {x 3 i[2(l - e 2 ™ /3 ) + (1 - e~ 2 ™ /3 )] 
+ x 32 [2(l - e - 2 W3) + (1 _ e 2*i/ayy 5z + 3^ = . 

(26) 

To solve this one can either use the quadratic formula, or 
note that if we neglect the Sz term this gives Sz ~ 1/ V^V, 
and then making use of this we see that both the constant 
term in the equation and (Sz) 2 are of order 1/N, while 
the term linear in Sz is of order I/./V 3 / 2 , i.e. smaller than 
the other two, so that it represents an even smaller order 
correction. This implies that to lowest order we have 



Define 6q so that 



'2i6 a 



(27) 



(28) 



which implies that z = — e ±2l6 ° and that 28q — \/3x~3~. 
The eigenvalues of U corresponding to these values of z 
are ±ie ±ie °. 

We can now find the cigcnstates corresponding to these 
four eigenvalues by making use of the relations 



a ■, 



V 



A 2 e"^ + 1 



(29) 



where A is one of the eigenvalues, and 77 is a normalization 
constant. This gives us the four eigenvectors: for A = 

ie l8 ° 



z 3 + z 2 



l = (z 2 -l)(z + I)=0, 



(23) 



which has a double root at z = — 1. Let us now calculate 
corrections to this root. First we set z = — 1 + Sz and 
substitute it back into the equation. We then keep terms 
of up to second order in small quantities, with both Sz 
and anything of order n 3 /N being considered small. To- 
ward this end, we set £31 = n 3 i/N and X32 = n 32 /N, 
and note that 
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:+s 31 (l_e*rf/3) +a . ffl(1 _ e -a*i/3) i (24) 



^ \j,0)+ie i °o\0,j) 

2-^< l _ e 2i8 

i=i 

^ b-,0)+z e ^ e ~ 2 "/3| 0;j) 

2-~i l _ e 2i0 O g-27ri/3 

j="31 + l 
(JV/2)+K32 



E 

j=(iV/2)+l 
N 

E 

j=(N/2)+n 32 + l 



|j,0)+»e i9 °|0,j) 

1 _ e 2i8 a 

\j,0) +ie l8 °e 27ri / 3 \0 1 j) 



I _ e 2i8 0e 2iri/3 



, (30) 



■5 



for A = -ie l8a 



following four vectors 



16) = v 



"31 I • n \ 



e 2i0 o 



JV/2 

E 

j="31+l 



b'.o) 



|o,j) 

ie ie °e- 27ri / 3 |0,j) 



I _ e 2ie g-27ri/3 



(Af/2)+«32 , . . i0 , ., 



1 _ e 2ie 

|j, 0) -ie iflo e 27 ™/ 3 |0,j) 



j=(JV/2) + l 
JV 

E] i _ , 2/0, l( _>ti ; :s 

j=(AT/2)+n 3 2+l 



(31) 



for A = ie- l9 ° 



^|j,0)+ie-«°|0,j) 



e -2i0 o 

|j, 0) +ie- 4e °e- 27ri / 3 |0,j) 



_i=i 

W/2 

J="31+l 



li,o) 



(JV/2)+n 32 

E 

j=(JV/2) + l 
N 

E 

j=(AT/2)+n 32 +l 



and, finally, for A = — ie" 



"°|0,j> 



1 _ e -2iB a 

\j, 0) + ie -""( 



'" ' :!7ri ^ 3 |0, j) 



1 _ e -2i6»o e 27ri/3 



(32) 



| Mj,0)-» e-*°|0,j) 



g-2ie 

|j, 0) -ie- 4e °e- 27ri / 3 |0,j) 



.7=1 
W/2 

! _ r -_'/'0 l(r -_'7r/ 

J=n 3 i+1 



(JV/2)+n 32 

E 

j=(iV/2) + l 
AT 

E 

i=(7V/2)+« 32 + l 



|j,0)-ie- i9 °|0,j) 

1 _ e -2i8 



\j,0) -ie- w «e 27Ti 



/3 |0,j) 



I _ e ~2i8 0e 2iri/3 



(33) 



We find that for these vectors to be normalized, we must 
have j] = (1/2) ypV. 

Now in order to help make sense of this, define the 



^(|j,0)±z|0,j» 



(N/2)+n 32 

+ Yl (li,o>±i|o,j» 

j=(JV/2)+l 



(34) 



and 



icr ) > = »7v^ 



g |i,0)±ie- 2 -/3| , j ) 

J = «31+l 



iV 



E 

i=(AT/2)+« 32 + l 



^ _ e -27ri/3 

|j, 0)±ie 2 ^/ 3 |0,j) 



1 _ e 2iri/3 

Now we have, approximately, that 

ia) = ^(ic 1 (+) ) + id +) )) 

l6) = ^(-|C 1 (+) ) + ld +) >) 

i6) = ^(ic 1 ( - ) ) + id- ) )) 
ie 4 ) = ^(-ic 1 ( - ) > + icr ) >). 



(35) 



(36) 



Now let us consider for the initial state of the walk the 
state: 



ivw) = -/=Ei j, °}- 

Defining the state 



(37) 



(38) 



we find that | (Co I V'init) I = i.e. these states have a 

substantial overlap. Noting that 



ICo) = |][>} 



(39) 



fc=i 



we see that after m steps 

+(-ie ie °) m \^) + (-ie- ie °) m \U)] 
= -±= {i m [ism(me )\({ +) ) + cos(m0 o )\( { 2 +) ) 

+ {-i) m [isin(m0o)|Ci _) ) +cos(m0 o )|C^ ) )]} 



(40) 
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What this means is that when m#o = 7r/2, |Co) becomes 
a state with the particle located on the edges with <fr = 0. 
Because the overlap between t/ T "'|V'i™t) and U m \Qo) is the 
same as the overlap between \ipmit) and \Co)> if we mea- 
sure the location of the particle in the state U m \ipi n it), 
we will find that it is, with probability 3/4 on one of the 
edges connected to a vertex with = 0. Note that in 
order to find edges with = with this procedure, we 
do not have to know which edges have — 2ir/3 and 
which have — —2ir/3. Also note that unlike in the 
usual Grover algorithm, the probability to find a vertex 
with = does not go to one as N — » oo. 



III. ORACLE VERSION 

Now let us show that we can rephrase our walk search 
problem as an oracle search problem. We consider the 
situation in which the function that the oracle evaluates 
is multivalued. The function /(j) with j = 1, 2, • • • ,N 
can take on integer values between and (d — 1), where 
d is an integer greater than 2. It is at j = jo, while it 
takes on one of the integer values between 1 and (d — 1) 
at other values of j. We are interested in finding jo. We 
will do this in two different ways. The first reduces to a 
standard Grover search, and the probability of finding jo 
using this method is close to 1. The second results in a 
search that is analogous to what happens in the quantum 
walk on the star graph, and in that case the probability 
of finding jo is a decreasing function of d. 

Let us first show how the oracle for a multivalued 
function can be used to perform a standard Grover 
search jUll^]. Our oracle has the action 0\j)\\k)2 — 
|j)i|fc© /(j))2, where |j)i is the state of the input regis- 
ter that contains the argument of the function, |fc)2 is a 
qudit state, and ® denotes addition modulo d. We need 
to manufacture an operation that flips the sign of | jo) but 
leaves other states |j), where j ^ jo alone. The operator 
0- x (h ®h- 2/i ® |0) 2 (0|)O, where I m for m = 1,2 is 
the identity on space m, does exactly that. We have 

0~ x (h ® h - 2h <g> |0) 2 (0|)O|j) 1 |0) 2 

|j)l|0>2 j^jo 

-b">i|0) 2 3= jo ' 



(41) 



because jo is the only value of j for which /(j) is 0. Note 
that because O d is the identity, we have that O^ 1 = O d ^ 1 
so that we can realize O -1 simply by applying O d — 1 
times. If we now apply this operator followed by the usual 
inversion-about-the-mean operator, we have the standard 
Grover iteration. This will result in a probability of close 
to one of finding jo. 

That is not how our quantum walk worked, however. 
There, we had d different phase shifts and different values 
of /(j) led to different phase shifts. Let us now formulate 
an oracle problem in which that happens. 

We assume for simplicity that /(j) = at only one 
particular value of j, say jo- To search for jo, we 



adopt a st raig htforward extension of the standard Grover 
algorithm 21 , [22j ■ The starting point of our search is to 
take a d x d matrix 



/I 



1 

1 P 

1 f3 2 



1 

Q4 



\1 P d ~ l d -^ ■ 



1 \ 



(42) 



as a d-dimensional equivalent of the Hadamard operator 
H = H2 = { 1 t ) , where 



1 -1 



P = e l27I / d {/3 d = 1), 



and 



1 + + (3 2 + ■ ■ ■ + p 



d-l 



0. 



(43) 



(44) 



With n qudits in the main register (d n = N) and one 
ancilla qudit in the secondary register, the search is then 
conducted in the following order: 

1. Initial preparation. The n qudits are prepared in 
state 



Cd—l \ ® n N 



(45) 



and the ancilla qudit in state H d \l) = 

2. Oracle transformation. The n qudits along with the 
ancilla qudit arc directed to an oracle in which the 
transformation \ j)\k) — > |j)|/(j)®^) is performed 
, where |j) and \k) represent the state of the n 
qudits in the main register and of the ancilla qudit 
in the secondary register, respectively. When the 
ancilla qudit is in state X)fc=o fi k \k), the oracle 
transformation is represented by 



d-i ^ 



Vd 



fe=0 



d-l N 



(46) 

The oracle transformation can thus be written, in 
the iV-dimensional Hilbert space of the n-qudit 
state, asO = £jli/r /( %)01- 

3. Inversion about average. This is achieved by 



N N 



N 



(47) 



1=1 j=i 



3=1 



4. Repeat the operation G = DO until the state of 
the n qudits in the main register comes sufficiently 
close to the desired state |jo). 
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5. Measure the state of the n qudits in the main reg- 
ister to get the solution. Evaluate /(j) to check 
that the solution is correct. If the solution is cor- 
rect, then the search succeeds. If not, go back to 
(1) and start the search over again. This last step 
(5) is necessary because, as will be shown later, the 
probability for the search to produce the correct 
solution is less than 1 for d > 3. 

In the TV-dimensional Hilbert space, the N x N matrix 
G = DO is given by 



N N 



fti) 



i=i 4=1 



i>0'i-$>" /w b'M- 

4=1 



After applications of the operator G k times, the state of 
the n qudits in the main register becomes \ipk) = G k \ipo), 
where \ipo) — 77= E4L1 \J)- The probability to find the 
desired state \jo) after k iterations is then given by 



P k = \{jo\G k \iM\ = 



N 



3=1 



where Zj's and \xj)' s are eigenvalues and eigenvectors of 
the matrix G. The probability P k can thus be deter- 
mined by calculating Zj's and |xj)' s - When d = 2, the 
method described here reproduces the results of the stan- 
dard Grover search with a two- valued function. 

Let us now find the eigenvalues z k (k = 1,2, •• • , N) 
and corresponding eigenvectors \xk) of the N x N matrix 
G given by Eq. (|48|) . where 



(50) 



G\xk)=z k \ Xk ), fc = l,2,-.. ,N 



The eigenvectors \xk) can be expressed in a linear super- 
position of the basis vectors | j) as 



N 

\Xk) = ^2c k j\j) 

3=1 



(51) 



Substituting Eqs. (gSJ) and (J3TJ) into Eq. ((3UJ| and solving 
the resulting equation for Ckj , we obtain 



2 c 

W k 



where 



N 

4=1 



.7 = 1,V 



.JV 



-/(4) 



Cfe4 



(52) 



(53) 



Substituting our expression for Cjk back into the above 
equation, we obtain 



N 1 



_ AT 
/3/0)z fe " ~ 2" 



(54) 



This is identical to Eq. (6). We note that one application 
of the Grover operator G is equivalent to two steps of 
the walk, and therefore the eigenvalue z of G is equal 
to the square of the eigenvalue, A 2 , of the operator U . 
We proceed to solve Eq. (l54l) for the case where there is 
one j, i.e., jo, at which f(j) = 0. We further assume for 
simplicity that the integer values between 1 and (d — 1) 
of the function /(j) are distributed evenly among other 
(N—l) j's. There is thus one j, i.e., jo, for which f(j) = 
and f-i j's each at which /(j) = 1, 2, • • • , d - 1. We 

assume that -^3r is an integer. Note that if we choose 
N = d n , where n is a positive integer, then 



JV-l 



(48) integer. Using the relation 



1 + x + x 2 H h x^ 1 = 



l-x d 
1-x 



is an 



(55) 



with x = —ftizk and noting that ft — 1, Eq. (|54|) can 
be rewritten as 



1 



N - 1 



(49) l + z k \d-ljl- (-z k 



1 



d-l 

d-l-^(-Zfe) 

4=1 



_ TV 
(56) 



d- 1 



N-d 
N 



d-l 



+ 1 = 0. (57) 

i=i 



In the limit N ^ 00, the above equation becomes 

d-l 



T )E(-^) 3 + 1 = - 



(58) 



4=i 



This equation, because it has real coefficients, will have as 
its roots either real numbers or complex conjugate pairs. 
By examining the above expression and its derivative, 
both of which must be zero at a double root, one finds 
that the only double root of this equation is z k = — 1. 
This can be seen as follows. Setting the derivative of the 
above equation equal to zero and then multiplying by 
-z k gives 



(59) 



Substituting into this equation for (— Zk) d from Eq. 
we find 



d-l 



d- 



4=i 



J 



(~Zk) j 



(60) 



Noting that 



2 , 



4=1 



(61) 
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we see that the only way Eq. (|60|) can be satisfied is if 
all of the terms add in phase. The only way that can 
happen is if zu = — 1, and so, as stated above, that is the 
only double root. 

Going back to the equation for finite N , we find that 
this root splits into two roots. In order to obtain an 
expression for these two roots, we substitute Zh = — 1 + Sz 

into Eq. (|57| and find Sz ps ±i^Jj^-^^. We designate 

these two roots as z\ and zi. We can now write 



where 



Z\ ps — 1 + i \ 



12 



N(d+1) 



Z2 ps — 1 — i\ 



12 



N(d+1) 



-i28 



= -e z29 



sin(20) ps 29 



12 



N(d + 1) 



(62) 
(63) 

(64) 



The finite N corrections of these roots are of order N^ 1 / 2 , 
whereas the finite N corrections to the other roots of Eq. 
(l58l) are of order AT -1 . 



Let us now find the eigenvectors |xi) and \xn) corre- 
sponding to z\ and z 2 . Using Eqs. (jST) and (f5"2"j) . we 
write 



Vi 



N 

E 3-/0') 



-i26 



\J) 



(65) 



and similarly for |xa}- The normalization condition 1 = 
(Xilxi) determines the constant 771. Using again Eq. (|55|) 
with a; = /3~ J e , the normalization condition is reduced 
to 



2(1 -cos 26>) V d - 1 / 2[1 - 



N - 1 



1 



cos(2rf6»)] 
I 



d-l 



[d(d-l)-2j2(d-j)cos(2j6)] 

3=1 



(66) 



,2 6 



For (n/N) <C 1, we have that cos2n6> ss 1 — n wpjji 
we obtain 771 ps ^/ N ^ +1 ^ , and thus 



IXi) 



and similarly 



1X2) 



6 



_ W 1 



iV(d+ 1)^/3 



/ JV(d+ 1)^/3 
Other eigenstates are given by 



_ N 1 

T)E g-/0')- e ^ lj) ( 68 ) 



Xi- 



>=^E fl _ /(j L fc = 3,4,.-.,JV (69) 



and the normalizations for these states are 

1 2 



JV 



/3-/Cj) + Zfe 



1. 



(70) 



For any j and fc, we have |/3 ^0') -f Zk\ < 2. Therefore, 
?7fe's should satisfy 



9 iV 
I%| 2 T <1 



(71) 



or 



m < 



N 



(72) 



The probability Pk to find the desired state after 
fc iterations is given by Eq. (l49l) . In order to evaluate 
this, we begin by noting that 



O'olXi) 



-i^=, O0IX2) 



1 



'V2' 



(73) 



which implies that (j'olXfc) ~ 0- Therefore, we have that 

Pk « I^O'olxiXxil^o) +40o|X2)(X2^o>r (74) 
From Eq. (|67[) . we have 

_ W 1 



(Xi|V>o) 



6 



AT 2 (d+ 1) ^ /3/0) - e' 



(75) 



and similarly for (x2|^o)- Evaluating the summation in 
in this equation using Eq. (|54[) . we obtain 



(xi|^o) = (x2|^o)«-y^ TT y 
Making use of these results, we finally obtain for Pk 



(76) 



P k 



d+1 



sm 2 (2k0). 



(77) 
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The maximum probability is given by P ma x ~ and 
the number of iterations to reach the maximum proba- 
bility is k max ps « f y^ ^^a" 1 " ■ Since the maximum 
probability decreases as c? increases, it gets increasingly 
important to check in the end that the solution produced 
by the search is indeed correct. This can be done cither 
classically or quantum mechanically. When d = 3, we 
have Pmax ~ f , in agreement with the result presented 
in the previous section. Although the maximum prob- 
ability decreases with d, the number k max of iterations 
required to reach the maximum probability still scales as 
y/~N for d > 3. This indicates that the quantum walk 
search still provides quadratic speedup, the best one can 
gain from a quantum search algorithm |23l - l25| . for d > 3. 

So far, we have assumed that there is a single match. 
We now briefly comment on the situation where there 
are multiple matches [24| . If the number M of matches 
is small compared with the total number N, only a triv- 
ial modification is needed. Since the algebra involved 
is straightforward, we here state the result only. The 
maximum probability P max remains the same as long as 
M N, but the number k max of steps to reach the 
maximum probability decreases with M. In fact, we ob- 
tain Pk « sin 2 (2k9) , where sin 9 = N ^+i) > ano - 

therefore k max ~ \J ^^^ - The number k max scales 

as 1/yM. Classically, when there are M matches, one 
would expect to find a match M times faster than when 
there is a single match. One thus expects that classically 
the number k max to scale as 1/M. It is then natural to 
expect that this number would scale as 1/VM for the 
quantum search which is quadratically faster than any 
classical search. If the number M of matches is not small 
compared with N, a different search strategy from that 
described so far should be employed because a match can 
be found with a reasonably high probability after only a 
small number of iterations. In general, the state of the 
qudits in the main register changes rapidly, and it seems 
difficult to find a general rule here. One simply needs to 
find the best strategy for a given situation. 

IV. CONCLUSION 

We have studied a quantum walk search on a star graph 
and an analogous oracle search problem. The quantum 



walk is trying to find the set of vertices that reflect the 
particle with a particular phase, and the oracle search is 
trying to find the input values for a multivalued function 
that produce a particular output value. We have shown 
that the maximum probability P ma x for the search on a 
star graph to find the correct vertex decreases with d, 
the number of different phases with which the particle is 
reflected from outer vertices. We have also shown that 
the number k max of steps needed to reach the maximum 
probability increases with d. These two effects combine 
to increase the average number of steps required to be 
taken in order to find the correct vertex. Despite that, 
however, this number still scales as yN and the quan- 
tum walk search still gives the quadratic speedup over a 
classical search for any value of d as long as N S> d. 

We examined two different oracle searches. One mim- 
icked the quantum walk, and for this search the prob- 
ability of finding the inputs that produced a particular 
output decreased with d, just as for the quantum walk. 
However, we also showed that the oracle can be used in 
a different way so that the probability of finding the de- 
sired elements is close to one. This indicates that the use 
of the oracle is more powerful than the quantum walk 
search we studied. 

Finally, we would like to note that it has been possi- 
ble in some cases to improve the maximum probability 
of finding a marked element in a quantum walks search. 
The original quantum- walk search on the hypercube due 
to Shenvi, Kempe and Whaley found the marked ver- 
tex with a probability of less than 1/2 [26]. Recently, 
methods have been found to improve this probability [27| ■ 
Whether similar methods could be used to improve the 
search on the star graph is a subject for further research. 
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